Impact of non-standard neutrino interactions on future 
oscillation experiments 
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Abstract We study the performance of reactor and superbeam neutrino experiments in the pres- 
ence of non-standard interactions (NSI). We find that for some non-standard terms, reactor and 
superbeam experiments would yield conflicting result in the &13 determination, while in other cases, 
they may agree well with each other, but the resulting value for S13 could be far from the true 
value. Throughout our discussion, we pay special attention to the impact of the complex phases of 
the NSI parameters and to the observations at the near detector. 
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1 Introduction 

In lepton flavour physics, one of the most anticipated 
upcoming results is the measurement or constraint of 
the mixing angle ^13 by reactor and accelerator neu- 
trino experiments. Within the context of standard os- 
cillations, the physics goals of these types of experi- 
ments overlap partly, but if physics beyond the stan- 
dard model should affect neutrino oscillations, we may 
obtain independent information from them because 
they observe the different oscillation channels, and new 
physics can therefore affect them in differing ways. 
In this talk, we will study the question of how non- 
standard effects could modify the results of reactor and 
superbeam experiments, and will specifically address 
the following questions: 

— What will it mean if reactor and accelerator exper- 
iments yield conflicting results? 

— If they give consistent results, how rigid is the stan- 
dard oscillation interpretation? 

In order to approach this subject in a model indepen- 
dent way, we will introduce effective four Fermi in- 
teractions (usually called non-standard interactions or 
NSI [T]|2] ) , which are assumed to be induced by physics 
beyond the standard model at some high energy scale. 

For NSI mediated by charged currents, the neu- 
trino states involved in the production and detection 
processes will not be pure flavour states, but can be 
described as mixed-flavour states, which we call ji/") 
for the source state and (v'^l for the detection state. 
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They are defined as [T] 



7— e,/j,,r 



(1) 

(2) 



Here, the index a stands for the flavour of the accom- 
panying charged lepton in the respective process. 

We can also introduce NSI mediated by neutral 
currents, which would create an additional, possibly 
flavour-violating matter potential of the form [5] 



(V^NSl)/3a = V2GFNee 



Pa- 



(3) 



The e parameters represent the ratio between the 
NSI amplitude and the standard amplitude. For exam- 
ple, e^^ and ejlg, are the degrees of "contamination" of 
the source and detection states with "wrong" flavours, 
and e^Q, denotes the ratio between the non-standard 
and standard matter effects in the propagation Hamil- 
tonian. Taking into account the NSI shown in Eqs. ([T]) 
- ([3]), we can calculate the oscillation probability for 
vf,, vi as 



f/3 



(4) 



Here, i?so is the standard neutrino propagation Hamil- 
tonian, which is parameterized by the standard oscil- 
lation parameters. Although the Lorentz structure of 
the effective NSI operators is not fixed by experimen- 
tal constraints, we will only consider {V — A){V — A) 
type interactions in this talkQ Under this assumption. 



^ The case of NSIs with different Lorentz structures is 
discussed in Ref. |3]. 
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Figure 1. Examples for a discrepancy (left) and a common offset (right) in the standard oscillation fits for Double 
Chooz and T2K. In the left panel, corresponding to the NSI parameter , the best-fit value of Double Chooz (vertical 



solid line) is located far from the best-fit value of T2K (coloured diamonds). In the right panel, with 
fits coincide with each other but have a common offset from the true value (denoted by the black star). 
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we can treat the e parameters as energy- independent. 
Moreover, the following relation between the NSI in 
the neutrino beam source and those in the detection 
process holds: 



(5) 



We will use this relation as a constraint in the numer- 
ical simulations which we are going to present in the 
next section. 



2 NSI-induced offsets and discrepancies in 
^13 fits 

To fully assess the high-level consequences of NSI for 
realistic reactor and superbeam experiments, we have 
performed numerical simulations using the GLoBES 
software [4,. In this talk, we will concentrate on only 
two experiments, namely T2K and Double Chooz. To 
analyse the simulated data, we define the following 
function 



channel bin 

=min J2 H 



Priors, 



|7V„■(A'™^e'™^)-7V,,(A,e = 0)|2 



iVij(A'™,e*™) 



(6) 



where Nij denotes the number of neutrino events in 
the z-th energy bin for oscillation channel j, the vector 
A = (^12, 6'i3, 6123, Sep, Am2i, Amli,b) contains the six 
standard oscillation parameters and the systematical 
biases b, and e represents the non-standard parame- 
ters. For the following plots, has been marginalised 
over all standard oscillation parameters (except, of 
course, those for which we are interested in confidence 
regions rather than best-fit points) and over all sys- 
tematical biases. We keep the NSI parameters fixed 
at in the fit because we want to study the impact 



of NSI on a standard oscillation fit. Unless indicated 
otherwise, we calculate the simulated event rates us- 
ing the following reference values (true values) for the 
standard oscillation parameters [9] , 



sin^ 26%'''' = 0.84, 



sin^ 261*™'= 



sin^ 26**™" = 1.0, 

ctrue 
Or 



0.05, d}^^" = 0, 

(Z\m^i )"■"'= = 7.9-10^5 eV^ 



= 2.5 • 10"^ eV^ 



(7) 



We assume the true hierarchy to be normal. Although 
we keep S^p^ fixed, this does not limit the general- 
ity of our results, because the oscillation probabilities 
typically depend only on combinations of Sqp and the 
phases of the NSI parameters. Thus, a variation of Sqp 
has the same effect as a variation of arg(e). 

In our simulation of T2K, the Super-Kamiokande far 
detector and a 1.0 kton water Cerenkov near detector 
are simulated separately. We introduce a common 10% 
uncertainty on the neutrino flux and a common 20% 
error on the number of background events in the Ve ap- 
pearance channel. Also for Double Chooz, we simulate 
the near and far detectors separately. We introduce a 
2.8% correlated flux normalization error as well as un- 
correlated 0.6% fiducial mass errors, and a bin-to-bin 
uncertainty of 0.5%. More details on the input param- 
eters of our simulations are given in Ref. [3] . 

Two exemplary outcomes of the simulation are plot- 
ted in Fig.[T] In both panels, we have chosen the "true" 
parameter values sin^ 20i3 = 0.05 and Sep = tt, as 
indicated by the black star. The vertical black line 
and the coloured diamonds represent the best-fit val- 
ues for Double Chooz and T2K , respectively, while 
the shaded vertical band and the coloured contours 
are the corresponding 90% confidence regions. Blue 
symbols and lines stand for the T2K normal hierarchy 
fit, and dashed magenta ones are for an inverted hi- 
erarchy fit. The NSI contributions are assumed to be 
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Figure 2. Distribution of tlie best-fit values for d\z in T2K and Double Chooz in the presence of the NSI parameters 
e^^ and i^^na)* ■ The black stars indicate the assumed true values of sin^ IQn which are taken to be 0.05 and 0.01. 



= 0.5e-''^/2 in the left hand plot and = 4e = 
0.05 in the right hand plot. In the first case, the best- 
fit value of Double Chooz conincides with the "true" 
value because reactor experiments are sensitive nei- 
ther to standard nor to non-standard matter effects. 
On the other hand, the ^13 fit of T2K deviates signif- 
icantly, indicating that e™ has a large impact on this 
experiment, even though standard matter effects are 
only sub-dominant due to the relatively short base- 
line0 The discrepancy is so large, that the best-fit 
value of T2K is almost ruled out at 90% confidence 
level by the reactor measurement. In the case shown 
on the right hand side of Fig. [H both experiments 
agree well with each other, but their fits suffer from 
a common NSTinduced offset, which even leads to an 
erroneous exclusion of the true 6*13. 

Let us now proceed to a more systematic analy- 
sis, which we present in Figs. [2] and [3l In these fig- 
ures, we plot the distribution of the best-fit ^13 val- 
ues of both experiments for different values of the NSI 
parameters. Each curve corresponds to fixed |e|, and 
varying arg(e), with dark red curves denoting |e| 0, 
and yellow curves corresponding to |e| equal to its cur- 
rent upper bound, which is 0.1 for e^^ [6_ and 0.7 for 
eS- 0- If the function for one of the experiments 
exceeds 3a, the colour of the curves has been changed 
to grey, regardless of |e|. This indicates that in these 



^ Standard and non-standard matter effects become 
more significant at baselines of 0(10'^) km, which are, for 
example, considered for a neutrino factory. For recent stud- 
ies of NSI in this context, see e.g. Refs. [71|8]. 
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Figure 3. Distribution of the best-fit values for &13 in T2K 
and Double Chooz in the presence of the NSI parameter e^. 



cases, a Sct discovery of the respective non-standard 
effect is actually possible. The black stars indicate the 
chosen "true" ^13. Let us briefiy comment on the dif- 
ferent cases shown in Figs. [5] and O If e^g and ef^ are 
introduced (upper left panel of Fig. [2]), they modify 
the De flux in Double Chooz, and lead to a wrong mea- 
surement of the Ve contamination in the T2K beam 
by the near detector. Therefore, T2K is affected indi- 
rectly. The correct treatment of the near detectors is 
crucial here to fully assess the impact of NSI. Similarly, 
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Figure 4. Scatter plot of the Oi^ best-fit values of T2K 
and Double Chooz in the presence of random combina- 
tions of NSI parameters. It is obvious that NSI can in- 
duce ostensible discrepancies between the two experiments 
(off-diagonal points), or common offsets (close-to-diagonal 
points) which would lead to consistent, but wrong results. 



the effect of e^^ and e^^ can only be fully understood 
if the near detectors are taken into account [3 . The 
upper right panel of Fig. [2] shows that e^^ and ef^ are 
particularly dangerous because they lead to a common 
offset in the reactor and in the superbeam (the best 
fit values are distributed along the diagonal) . We have 
already seen an example of this effect in Fig. [TJ The 
NSI parameters e^^, and e^^ affect only T2K, hence 
the distribution of the corresponding 6'i3 fits is con- 
strained to horizontal lines in our plots. Among the 
non-standard matter elfects, only can give a sig- 
nificant effect, because the other parameters are ei- 
ther subdominant in the oscillation probability, or are 
already well-constrained experimentally. The effect of 
e™ is a possibly large shift of the 6*13 fit in T2K, as 
shown in Fig. [3l 

So far, we have only considered situations in which 
one NSI parameter is dominant, and all the others are 
negligible. In contrast to this. Fig. 2] shows the dis- 
tribution of the fits in situations where many NSI 
parameters are present simultaneously. For each point, 
the moduli |e*'^'™| were chosen randomly on a logarith- 
mic scale between 10~® and the respective experimen- 
tal upper bounds [5l[6], while the phases are linearly 
distributed between and 2tt. The colour coding in- 
dicates the value for each point. We can see from 
the plot that, for part of the parameter space, the NSI 
effect can actually be discovered, i.e. the standard os- 
cillation interpretation does not give a good fit. How- 
ever, there are also a lot of points with a low x^i and 
among these are many which exhibit a clear discrep- 
ancy between the 6'i3 fits of Double Chooz and T2K, 
and others which correspond to a common offset of the 
fit values. 



3 Conclusion 

In this talk, we have discussed the impact of non- 
standard neutrino interactions on reactor and accel- 
erator experiments, in particular on Double Chooz and 
T2K. We have introduced flavour violating four-fermion 
interactions with {V — A){V — A) type Lorentz struc- 
ture, and considered the impact of these NSI oper- 
ators on the neutrino production, propagation, and 
detection stages. We have performed detailed numer- 
ical simulations, taking into account parameter corre- 
lations, degeneracies, systematic errors, and a realistic 
treatment of the near detectors. We have found that 
NSI can have a sizeable impact on future experiments 
if the corresponding couplings are close to their cur- 
rent upper bounds, and if the complex phases do not 
conspire to cancel them. There are scenarios in which 
a clear discrepancy between Double Chooz and T2K 
arises, but we have also found situations in which both 
experiments would agree very well, but the derived ^13 
has a significant offset from the true value. To detect 
this kind of problems, a third experiment, complemen- 
tary to the other two, would be required. Our discus- 
sion shows that reactor and superbeam measurements, 
which might seem to be redundant in the standard os- 
cillation framework, turn out to be highly complemen- 
tary once non-standard effects are considered. 

More details on the impact of non-standard inter- 
actions on reactor and superbeam experiments can be 
found in Ref . [3J , in which we also discuss the discovery 
potential for NSI, paying specific attention to its de- 
pendence on the complex phases. Moreover, Ref. 3 pro- 
vides an intuitive understanding of the impact of the 
different e parameters, and contains analytic formulae 
to substantiate the numerical results. 
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